Theory of anisotropic Rashba splitting of surface states 
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We investigate the surface Rashba effect for a surface of reduced in-plane symmetry. Formulat- 
ing a kp perturbation theory, we show that the Rashba splitting is anisotropic, in agreement with 
symmetry-based considerations. We show that the anisotropic Rashba splitting is due to the admix- 
ture of bulk states of different symmetry to the surface state, and it cannot be explained within the 
standard theoretical picture supposing just a normal-to-surface variation of the crystal potential. 
Performing relativistic ab initio calculations we find a remarkably large Rashba anisotropy for an 
unreconstructed Au(110) surface that is in the experimentally accessible range. 

PACS numbers: 71.15.Rf 73.20.At 75.70. Tj 



i 

C 

o 
o 



> 
o 

O 

o 



X 



Metallic surfaces often exhibit Shockley-type surface 
states located in a relative band gap of the bulk band 
structure, and forming a two-dimensional electron gas. 
One of the most intriguing manifestation of spin-orbit 
coupling (SOC) at surfaces is the splitting of these surface 
states, known as Rashba splitting P, 0|. Such Rashba 
splitting was observed via photoemission by LaShell et 
al. for the L-gap surface state at Au(lll) and ex- 
plained theoretically in terms of a tight-binding model Q 
and ab initio electronic structure calculations Q , but 
several studies of the Rashba splitting were published in 
recent years on Bi(lll) and Bi/Ag(lll) [7j-|9|, as well as 
on Bi K Pbi_3;/Ag(lll), where atomic Bi p-orbitals lead 
to a more pronounced spin-orbit splitting 

Describing and controlling the Rashba splitting of sur- 
face states is crucial for spintronics applications. The 
famous Datta-Das transistor relies on the electric tuning 
of the Rashba splitting [l4| and the Rashba splitting is 
responsible for the spin Hall effect in two dimensions [15[ 
and the anomalous Hall effect [HI as well. 

The simplest way to understand the origin of the 
Rashba effect is to take nearly free electrons, confined 
by a crystal potential, V(r) = V(z), and having a plane- 
wave- like wave function, ipsM (r) = e lkr </> (z) \s, with \s 
some spinor eigenfunctions, and k the momentum paral- 
lel to the surface. The crystal potential V(z) obviously 
produces an electric field, E, perpendicular to the sur- 
face, which, in the presence of spin-orbit interaction leads 
to the following spin-orbit term in the effective Hamilto- 
nian, 



H R (k) = a R (k x a y - k y a x ) 



(1) 



called Rashba-Hamiltonian. 
Pauli matrices and a R = 



In Eq. ([T]), cr, denote the 
^/d 3 r|0(z)| 2 a z nr) is 



the so-called Rashba parameter. The eigenvalue prob- 
lem can then easily be solved, resulting in a splitting 
of the spin-degeneracy of the surface states, e± (k) = 
T^--k 2 ±af? Ikl , with m* the effective mass of the surface 
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electrons 0, ■ Clearly, the above dispersion is isotropic 
in k-space, hence we term it as isotropic Rashba splitting. 

Although real systems cannot be described in terms of 
free electrons, and for quantitative estimates of an the 
atomic structure of the potential needs be taken into ac- 
count [|| , the structure of the Rashba interaction, Eq. ([1} , 
is very robust for surfaces of high point-group symmetry 
such as CVu or C± v [lit ]. 

The situation is, however, quite different for surfaces 
(or points in the surface Brillouin zone) of reduced sym- 
metry. Such Shockley-type surface states emerge, e.g., 
around the Y point of the Surface Brillouin Zone of un- 
reconstructed and (2x1) reconstructed Au(110) surfaces, 
as revealed by recent high-resolution photoelectron spec- 
troscopy experiments [17] . In this case, the Civ point- 
group symmetry of the system not only implies the asym- 
metry of the effective mass, m* ^ m* (for the crystal axes 
see FigQ} but, in leading order in k, representation the- 
ory also predicts the following simple form of the effective 
Hamiltonian I19II. 



H(k) =e - 



h 2 k 2 x 
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The above expression can easily be justified by simple 
symmetry analysis, just by noticing that a v and —o x 
transform as p x and p y under the operations of the dou- 
ble groups of Civ and C^ v . From this observation it 
also follows that in case of Ca v point-group symmetry 
ctR.x = ctR.y must be satisfied, and the Hamiltonian ([T]) 
is recovered. 

Although the above form of the Rashba interaction has 
been predicted in Ref. (l9| . no microscopic theory has 
been constructed so far to support it. While previous 
ab initio calculations [I?], 18 1 did find a Rashba split- 



ting of the Au(110) surface state, they focused only on 
the dispersion along the TY direction, and therefore the 
anisotropy of the Rashba term remained unnoticed. In 
the present paper, we provide such a microscopic analysis 
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FIG. 1: (Color online) Left: Sketch of the fcc(llO) Surface 
Brillouin Zone. The dark area denotes the projection of the 
L-gap of bulk Au. Right: Structure of the surface energy 
spectrum in the absence of SO interaction, along the line k = 
(k x ,0). Surface states in the relative gap with k / can be 
built up from states indicated by the thick black lines and the 
black circle at k = 0. Note that k = corresponds to the Y 
point of the Brillouin zone,see Eq. @- 



for an Au(llO) surface with Civ point-group symmetry. 
First, constructing a k ■ p perturbation theory for the 
surface states we show that the above anisotropic Rashba 
structure appears naturally, and is due to the finite mo- 
mentum mixing of the bulk p states to the surface-state. 
We also perform ah initio calculations of the Rashba-split 
surface state of an unreconstructed Au(llO) surface and 
confirm with a high numerical accuracy that there is a 
large anisotropy in fc-space, a^ x ~ 5q^ 9 , in agreement 
with Eq. ([2]). The predicted anisotropic Rashba splittings 
turn out to be within the range of experimental accuracy. 

Bloch-states of Au(110) can be characterized by a sur- 
face momentum, and can thus be written as 



^q+k (r) = e tkr Q , k (r) 



(3) 



with the momentum k measured with respect to the mo- 
mentum Q associated with the Y point of the Surface 
Brillouin Zone. Here the functions 0Q,k (r) are lattice- 
antipcriodic in the x direction, while they are lattice pe- 
riodic in the y direction of the (110) plane, see Fig. 1. 
For any given momentum, k, there exist an infinite num- 
ber (continuum) of eigenstates, the energy of which (e k ) 
is determined by the condition that the states V'Q+k be 
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eigenstates of the Hamiltonian, H = ^— + V(r) + Hso, 
with Hso denoting the spin-orbit coupling, 



Hso (r) = 



(VV(r) xp)cr 



(4) 



As a consequence, the functions 0Q,k must satisfy the 
equation, 



P 

2m 



V(r) 



h 2 k 2 
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k • p + Hso (k, r] 



£k0Q,k(r) , 



:,k(r) 
(5) 



with Hso (k, r) being the effective SO coupling, 

n 



H so (k,r)=H so (r) 



Am 2 c 



( W (r) x hk) a . (6) 



Similar to Bloch wave functions, for any fixed momen- 
tum, k, (and for any value of Hso) the functions 4>QM 
form a complete set for functions having the previously- 
mentioned periodicity property. In the spirit of k • p 
perturbation theory, we can thus take the complete set 
of k = and Hso = solutions, satisfying 



2m 



4>i,ni (r) 



,^,n, (r) 



(7) 



and expand (j)Q,k in terms of these. Here we classified 
the solutions according to the four one-dimensional irre- 
ducible representations of the Ci v symmetry associated 
with the point Y, j £ {1, x, y, ccy}, and labeled solutions 
of a given symmetry by rij. As shown in Fig. [TJ the spec- 
trum contains a discrete surface state of s-symmetry and 
the projected bulk continuum forming the gap. Let us 
denote the k = surface state by <f> , and its cigenenergy 
by £q. Then states with k ^ but with Hso = can 
be expressed in terms of the states 4>i,ni by performing 
second order perturbation theory in k, which amounts in 
a surface state 
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with approximate dispersion 
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(i = x,y) (10) 



The index in e£ and |(/) k ) is meant to remind us to the 
absence of SO interaction. 

To obtain the surface states, |</> k ), we then carry out 
first-order perturbation theory with the SOC operator, 
Hso 5 using the states as a starting point. Keeping 
just contributions linear in k we get two terms to the ef- 
fective Rashba Hamiltonian. The second term in Eq. ([6]) 
gives rise to the usual isotropic Rashba model, 



TTlSO 



with an 



(k) = a R (e z x k) ■ a , (11) 
The term Hso m 



Eq. ([6]), however, gives also a finite contribution due to 
the admixture of p XjV states from the continuum and, in 
fact, this is precisely the term that leads to an anisotropic 
Rashba coupling, 



Tj ams 



(k) = (12) 

{4>i,ni \Pi\4>o) ((0o| a l^.^rti) &) + h.C. 
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where we defined the (axial)vector operator related to 
SOC, a = (W (r) x p). Using the symmetry of 

the unperturbed wave functions, a particularly simple 
form of the above anisotropic Rashba Hamiltonian can 
be obtained, 

H% nls (k) = \ x k x a y + X y k y a x , (13) 
with the coefficient A^ expressed as 

A = 2 Re((^,»JPx|^o) (</>o|« y \<fix,nj) , 14 » 

X m ^ £ - £ x ,n x 

and A y given by a similar expression. The structure of 
this term is identical to the one obtained by symmetry 
analysis. Decreasing the in-plane asymmetry of the po- 
tential (e.g., considering C^ v point group symmetry), the 
relationship, X x = —X y , is satisfied, i.e., the Rashba split- 
ting becomes isotropic. Nevertheless, even in this case, 
the anisotropic Rashba term, Eq. (TlB"]) . arising from the 
mixing of the surface state with bulk states, also con- 
tributes to the Rashba splitting. Thus we find that, up to 
second order in k, surface states of Au(llO) are described 
by Eq. ((2]), with un. x = (Xr + X x and OiR, v = -an + X v . 

To obtain a quantitative estimate of the parame- 
ters a R _ x / y and the induced Rashba splittings, we per- 
formed calculations of the surface states of unrecon- 
structed Au(llO) surface near the Y point of the Surface 
Brillouin zone, using the relativistic Screened Korringa- 
Kohn-Rostoker (KKR) method. Details of this method 
are described in Refs. 0, H^]. The computed disper- 
sion relations along the TY (x) and the YS (y) direc- 
tions are plotted in Fig. [2] The maximum binding en- 
ergy, eo ~ 370 meV, is by about 200 meV less than the 
measured value [13] and other theoretical values [13, EH ■ 
This deviation is mostly caused by the atomic sphere ap- 
proximation (ASA) and the angular momentum cut-off, 
(-max — 2, which resulted in some error for the determi- 
nation of the Fermi level and the vacuum potential. 

The nearly free electron-like, parabolic shape of the 
dispersion as well as the Rashba splitting being remark- 
ably different along the two directions is obvious from 
Fig. [2] and a detailed analysis confirms this impression: 
The numerical results are very well fitted by the disper- 

sions e ± (k) = s + ff + £J ± ^ n] ; ,ki ■ , 
obtained by diagonalizing the approximate Hamiltonian, 
Eq. ([2]), with the fitting parameters, to* = 0.11 to, 
to* = 0.32 to, aR >x = 0.8 eV A, and a^ y = 0.17 
cV A. The obtained effective mass along the YS direc- 
tion is in satisfactory agreement with the measured value, 
m y = 0.25 m. [ItJ The effective mass along TY, m x is 
only about one third of m y , which is the consequence 
that the states at the lower bulk band edge are mainly 
of p z and p x character (see Eq. ((9])). Note that the en- 
ergy separation of the surface state at the Y point is 0.8 
and 3.4 eV with respect to the lower and upper bulk 
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FIG. 2: (Color online) Dispersion relations of the Au(110) 
surface states at the Y point (|fe| =0) along the Yr and the 
YS directions. Symbols refer to the calculated data, solid 
and dashed lines to the fitted curves for £-(k) and e+(k), 
respectively. 

band edges, respectively, implying a strong admixture of 
"electron" states from the continuum below the surface 
state. 

One of the most astonishing results of these numeri- 
cal calculations is the remarkably large anisotropy of the 
Rashba parameters, an. x ~ 5an. y . In view of Eqs. (|14l) . 
this observation can also be explained with the absence 
of p y states at the lower bulk band edge. This result 
also correlates with the results of the effective mass: the 
smaller value of m x indicates a stronger admixture of p x 
states, also responsible for the stronger renormalization 
of a x . We remark that a^ x is even larger than an we 
calculated for the L-gap state of Au(lll), 0.57 eV A. 
This latter value is though considerably larger than the 
experimental one, 0.4 eV A (23|, which correlates with 
the theoretically computed effective mass, m* ~ 0.19 m e , 
being too small as compared to the experimentally ob- 
served value, to* ~ 0.25 m e . The computed Fermi wave 
numbers, kp= 0.160 and 0.189 A -1 , on the other hand, 
are almost in perfect agreement with the measured val- 
ues [23[. Nevertheless, based upon the discrepancy re- 
garding the value of the effective masses, we expect that 
our theoretical calculations for Au(110) somewhat over- 
estimate the Rashba parameters, ctR tX / y . 

The anisotropic Rashba coupling together with the 
anisotropic effective mass gives rise to two Rashba-split 
Fermi surfaces for the surface states, as shown in Fig. [3] 
The Rashba splitting along YS, Ak y ~ 0.017 A -1 is in 
the order of the experimental resolution (0.01 A -1 ) [l7T |. 
but the value Ak x ~ 0.026 A -1 implies that the Rashba 
splitting should be detectable experimentally along the 
TY direction. 

In Fig. |4] we also show the polar plot of the energy split- 
ting, Ae(k) = e+(k) — £-(k), between the two bands at 
the inner Fermi surface (e+(k) = sp). As a compari- 
son, the same quantity is displayed for the surface state 
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FIG. 3: (Color online) Rashba splitting of the Au(110) surface 
state at the Fermi level. Solid and dashed lines refer to the 
bands, £+(k) and £-(k), respectively. 
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FIG. 4: (Color online) Energy differences, Ae(k) = £+(k) — 
e_(k), for the Rashba-split surface state of Au(110) (solid 
line) and Au(lll) (dashed line) as a function of the polar 
angle <p = arctan -r 2 - , shown in units of degree around the 
graph. The magnitude of k was fixed to satisfy e+(k) = ef- 
The energy scale is indicated by the axis on the left. 

of Au(lll). Supporting the above implication, Ae(k) 
along YT for Au(110) is almost as large as the (isotropic) 
energy splitting in case of Au(lll). Furthermore, the 
extremely strong anisotropy of Ae(k), that could be in- 
ferred from angle-resolved photoemission experiments, is 
a clear fingerprint of the anisotropic Rashba effect dis- 
cussed in this work. 

In summary, we constructed a k • p perturbation the- 
ory for surface states in the presence of SO coupling, and 
derived a generalized Rashba Hamiltonian for (nearly 
free) electrons on metal surfaces. We found that in 



case of C2v symmetry, the Rashba interaction gets an 
anisotropic part in first order of k, which for Au(110) is 
found to dominate over the additional, well-known sym- 
metric term. The anisotropic Rashba term appears due 
to the mixing of the surface state with the bulk states 
for finite momenta. Even for surfaces of higher sym- 
metry, this mechanism (i.e., the corresponding term in 
k • p perturbation theory) gives a large contribution to 
the isotropic part of the Rashba Hamiltonian. Based on 
fully rclativistic first-principles electronic structure calcu- 
lations, we also demonstrated that a strongly anisotropic 
Rashba coupling should be experimentally observable for 
Au(110) surfaces. 
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